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In principle a quantum system could be used to simulate another quantum system. The purpose 
of such a simulation would be to obtain information about problems which cannot be simulated with 
a classical computer due to the exponential increase of the Hilbert space with the size of the system 
and which cannot be measured or controlled in an actual experiment. The system will interact 
with the surrounding environment, with the other particles in the system and be implemented using 
imperfect controls making it subject to noise. It has been suggested that noise does not need to be 
controlled to the same extent as it must be for general quantum computing. However the effects of 
noise in quantum simulations and how to treat them are not completely understood. In this paper we 
study an existing quantum algorithm for the one-dimensional Fano- Anderson model to be simulated 
using a liquid-state NMR device. We calculate the evolution of different initial states in the original 
model, and then we add interacting spins to simulate a more realistic situation. We find that states 
which are entangled with their environment, and sometimes correlated but not necessarily entangled 
have an evolution which is described by maps which are not completely positive. We discuss the 
conditions for this to occur and also the implications. 

PACS numbers: 03.65.Yz,03.67.Ac 



I. INTRODUCTION 

Simulating quantum systems with quantum systems 
is one of the primary reasons there is a great deal of 
interest in building a quantum computing device. The 
difficulty of simulating quantum systems on a classical 
computer, mainly due to the exponential increase of the 
Hilbert space with system size, was Richard P. Feynman's 
motivation for proposing the idea that a quantum system 
might perform this task much more efficiently [1 ] . Lloyd 
showed later that some quantum systems could be ma- 
nipulated to represent the evolution of other quantum 
systems using only local interactions [2]. 

There are many problems of interest in quantum me- 
chanics which have no known analytical solution. Thus 
for a wide range of physical systems simulation is a valu- 
able tool for solving quantum mechanical problems. Clas- 
sical simulation of such systems can quickly become in- 
tractable as the number of particles increases. The re- 
sources that are required to perform such a task increase 
exponentially with the size of the system. For example, 
in order to represent the state of N 2-state particles a 2^ 
vector is required and for its evolution the unitary will be 
a.2 N x2 N matrix However, only N particles would 

be necessary to simulate such a system [21 0]. In this 
sense, a quantum simulator is conjectured to provide ex- 
ponential speedup over classical simulation [5]. But that 
is not the only advantage; other problems such as the 
sign problem from Quantum Monte Carlo algorithms for 
fermionic systems, or the exchange-correlation function- 
als in Density Functional Theory [SI [7] will not be present 
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in a quantum simulation. Therefore, many difficult prob- 
lems in particle physics, condensed matter systems, and 
chemistry, among others, could be tackled [51 lol IBH2T)]. 

Quantum simulations have received a great deal of re- 
cent attention since they are feasible without the need for 
a universal quantum computing device. The question of 
the universality of Hamiltonians has been addressed to a 
great extent [2TH5U] and algorithms have been developed 
to simulate specific systems [1 H [m E3 ESSO] • In addi- 
tion, experiments have been designed and implemented 
[IB ] I4TH47] . However, a great deal of work remains to be 
done. Currently available quantum simulating devices 
have relatively few controllable particles. They are, af- 
ter all, quantum systems that inevitably interact with 
the surrounding environment and therefore are subject 
to noise. Just as with quantum computing, this is an im- 
portant issue when it comes to scalability. It is therefore 
necessary to study how the interactions affect a quantum 
simulation. 

The purpose of the present work is to study effects 
of noise in an existing algorithm proposed for a quan- 
tum simulation and to take away from this example as 
much general understanding as we can. The primary 
noise considered is prior unknown correlations or entan- 
glement within the system and between the simulated 
system and the environment. We study the evolution of 
different initial states, including ideal ones and states in 
which errors are present due to mistakes in preparation 
and/or interactions with particles in the system and find 
the dynamical maps that represent the evolution. The al- 
gorithm we explore was proposed and developed by Ortiz 
et al. 6J to simulate the one-dimensional Fano- Anderson 
model. To examine various behaviors of the system with 
initial correlations, we first provide a background for the 
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quantum simulation in Section IA" which focuses on the 



different sources of noise that can affect the experiments. 
Section |IB| provides a brief review of open system quan- 
tum dynamics, and discusses dynamical maps and their 
main characteristics, including requirements for positiv- 
ity and complete positivity; the purpose is to use dy- 
namical maps to describe general errors in simulations. 
Section [TT] contains a brief explanation of the algorithm 
used, including the modifications we made to represent 
noise in the system. Finally, our results, given in Section 
WJ are divided in two parts: those states for which the 
Bloch vector only has a component along the z direc- 
tion, and those which have some small component along 
x and a main one along z. We will also discuss why this is 
important. These two last subsections are subsequently 
divided into simulations performed with no external noise 
and simulations with noise. For the purpose of compar- 
ison, the parameters of the system were obtained from 
Ref. El and were used for all the considered scenarios. 



A. Quantum simulations 

There are two classifications of quantum simulators. 
The Universal Quantum Simulator (UQS) [48] (also re- 
ferred to as Digital [49] ) . is a quantum computer repre- 
sented by the standard circuit model with the set of uni- 
versal gates that act on a collection of two-state systems 
[2211501 151]. The term universal implies that the quantum 
computer would be able to simulate any arbitrary quan- 
tum system [52] which implies universal quantum compu- 
tation is possible. However, a fully functioning quantum 
computer has not been built yet. So researchers have 
designed and implemented devices consisting of smaller 
and controllable quantum systems specifically intended 
for simulations. This is the other type of quantum sim- 
ulators, referred to as Specialized Quantum Simulators 
(SQS) jlHl [55] or analogue quantum simulators US] [SS] 
since they are not able to be used to simulate any quan- 
tum device or computation. Rather, they are able to sim- 
ulate a smaller, but interesting class of physical systems. 
Examples of such systems include: ultracold atoms, ion 
traps, quantum dots, atoms in optical lattices, coupled 
cavities, photons, electrons floating on He films and NMR 
devices gl HH HS1 HS1 E21 S21 H31 IHl HH E31 HI]- In the 
SQS, universality for all quantum systems is not required, 
thus many interesting advances and simple simulations 
have already been performed [El UM 133 Ell US [HJ- 



Just as is the case with any other quantum system, 
unwanted interactions with an environment can have 
a detrimental effect on the outcome. Error correction 
and/or prevention is usually required for accurate imple- 
mentation. However, inaccurate unitary transformations 
are also a source of noise and the evolution of the sys- 
tem under a specific Hamiltonian is the main problem of 
interest [3j[60]. 

All steps, preparation, evolution and measurement, 



can cause some error EJ [T7] as well as unwanted inter- 
actions with other particles in the simulator, etc. It was 
initially suggested that decoherence in quantum simula- 
tions may not need to be treated in the same strict sense 
as in quantum computation 2j since noise in the simu- 
lating system might be able to be identified with noise 
in the simulated system. The nature of the interactions 
of the simulator with the bath may not be the same as 
those of the system of interest and thus error prevention 
techniques of some sort will almost certainly be required. 
These include error-correcting codes (QECC) [5TH5o] . 
decoherence free subspaces/noiseless subsystems (DNS) 
[oTHT2"] (see also [73] [73] for reviews), and/or dynamical 
decoupling (DD) [531 I75H85] . However, error correction 
means an increase in resource requirements, and this can 
represent a problem with scalability [3] [U [551 [551 ISZ] as 
well as efficiency. In addition, there exist algorithms and 
observables which have an inherent robustness to errors 
[88] , but this is not the case for all systems and all errors. 

It is sometimes, in fact, possible for the errors to be 
treated quite differently in the simulation of quantum 
systems. For example, sometimes it is possible to model 
some of the interactions of an open quantum system, as 
is done in Refs. [331 GUj, in which the bath is simulated 
as well. Also, Diir et al. propose an algorithm to gen- 
erate many body interactions |89) to study the influence 
of noise in the simulation of many body systems. Fur- 
thermore, an experimental setup to study open systems 
is proposed in [18] . In these cases, the environment is in- 
cluded, but external interactions will still be present and 
will affect the final outcome in an undesirable manner. 

Our work examines interactions within the system, 
where errors in preparation of the initial state give rise to 
errors. Obviously the initial state is important because 
the outcome of the simulation depends on it. Also, when 
errors are caused by initial entanglement, dynamical de- 
coupling cannot remove those errors since these controls 
rely on local unitary transformations to eliminate Hamil- 
tonian interactions with a bath. Local unitary controls 
cannot change the entanglement between the system and 
the bath. 

Experimentally, it has been observed that two different 
state preparation methods may not yield the same result 
and can have a profound effect on the outcome [90]. We 
observe the characteristics of the dynamical map, (which 
will be described more in detail in the next section) that 
describe the evolution of different initial states and de- 
termine their positivity or complete positivity. Until re- 
cently, discussions of the evolution of an open quantum 
system were limited to completely positive maps. How- 
ever, work by Pechukas 911 and more recently by Shaji 
and Sudarshan [92] have provided demonstrations that 
a map does not need be completely positive for the end 
result to represent a physical state. It fact, the map does 
not even need to be positive; it must only be positive on 
a given domain in order to possibly represent a physi- 
cal mapping. In certain circumstances dynamical maps 
can provide information about correlations in the initial 
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state of the system, which could provide useful informa- 
tion about the effects of noise and interactions in quan- 
tum simulations. Furthermore, there are many sets of 
operators in the operator-sum decomposition which give 
rise to the same map. This is true of completely positive 
maps [93l El] as wei l as maps which are not completely 
positive [9"5] . 



B. Noise in Quantum Systems, Completely and 
Non Completely Positive Maps 

The density matrix, or density operator, represents our 
knowledge of the quantum state of a system. In general 
any density operator must satisfy the following conditions 
in order to represent a physical state [96 : 



it is Hermitian, 



(1) 



p > 0, it is positive semi-definite, 

i.e. its eigenvalues are non-negative, (2) 

Tr(/o) = 1, it has trace 1, 

i.e. the sum of the probabilities is 1. (3) 

The evolution of a closed system is described by a unitary 
transformation, as 

m = [^(0), 

where U = exp iHt. It follows that 

pit) = Up(0)tf. 

The density operator is often written as an expansion 
of pure states 



= J2pj li) 01 > 



pure state. Thus any state on the surface of the Bloch 
sphere is a pure state. A mixed state is represented by 
a vector with | a |< 1. With this notation it is possible 
to have a visual representation of the quantum states at 
different times. 

A system S that is coupled to an environment E with 
Hilbert spaces T-Ls and He, respectively, can be consid- 
ered a larger isolated system whose initial state is de- 
scribed by psb(0). The time evolution of this system 
is then given by the joint evolution of the system and 
environment |92) 

p SE (t) = Up SE (0)Ul 

We are often only interested in the evolution of the sys- 
tem, S. Tracing out the environmental degrees of free- 
dom provides us with the reduced dynamics of the system 



p s (t) = Tr E [psE(t)} = Ti-e \Us E (t)psE(0)Ul E 



With the reduced dynamics of S, we can find the map 
that transforms the initial state p(0), into the final state 
pit). To obtain the "dynamical map" it is convenient to 
write the N x A density operator p as a TV 2 x 1 column 
vector that is transformed into another TV 2 x 1 column 
vector through the A 2 x N 2 supermatrix A 



Pr's'i*) = A r's'.rsPrs(0), 



(4) 



where A describes the most general evolution of p |104j . 
In matrix notation 



Ap. 



(5) 



Because p must be mapped to another positive p the 
following conditions are imposed on A |96j : 



A 



r s ,rs 



{A 



s r ,sr 



)* , is Hermitian, 



(6) 



where the pj are the probabilities associated to each of 
the states \j). If one of the probabilities is equal to 1 
and the rest are 0, then the state is pure. For two-state 
systems we can write the density operator in terms of the 
2x2 unit matrix and the Pauli operators, 



P : 



(1- 



where the coefficients ai are the projections along the x, 
y and z directions of the so-called Bloch vector. This 
provides a representation of the quantum state in the 
well-known Bloch sphere, which is a geometric represen- 
tation of the states of a qbits in terms of a sphere with 
radius 1. (For higher dimensional systems, this is referred 
to as the polarization vector, coherence vector, or gener- 
alized block vector. See [9"?lU03] and references therein.). 
The magnitude of the Bloch vector is constrained by the 

condition , /a 2 + a 2 + a 2 < 1, and | a |= 1 represents a 



x*x s A rs y s <y*,y s > > ,A is positive, (7) 



22 A rr y s ' — S r > s > ,A is Trace Preserving. (8) 

T" 

These conditions ensure the conditions Eqs. ((3]) on 
the density operator are satisfied. 

By interchanging indices of A, we obtain another TV 2 x 
A 2 supermatrix B [S3] 



B ' > = A , i. 

-^rr ,ss rs,r s 



(9) 



The 1 x A 2 rows of A become the A x N block matrices 
of B. The following conditions are imposed on B so that 
it represents a physical map: 

B rr > ss > = (B r , r , s T, B is Hermitian, (10) 
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x*y r >B rr ' ss 'X s y*, > 0, B is positive semi-definite, 

rsr' s' 

(11) 

/^rr'.rs' = Vs'i B is trace preserving. (12) 



As mentioned above, when the map is completely pos- 



itive the eigenvalues of B in Eq. ( 13 1 can be taken to all 
be positive. When they are, Eq. ( |13[ ) can be rewritten as 

p(t) = B [p(0)\ = ]T \ a ( aP (0)(l = £ C a p(0)Cl (14) 



From these we may write 



where C a — ^J\ a Ca- Eq. (14) is sometimes known as 



p(t) = B \p(0)] 



(13) 



If B is decomposed into its eigenvectors and eigenvalues, 
the action of the map can be represented as follows 

£[ P (0)]=]TA Q CaP(0)Ci, 



where A a £ R are the eigenvalues. The hermiticity of p' 
is guaranteed by the restriction given in Eq. ( 10 ) |104j . so 
that B must be Hermitian. The matrix A is required to 
transform p(0) into another Hermitian state p(t), but A is 
not necessarily Hermitian itself. The complete positivity 
of the map implies that the final state will be positive. 
The eigenvalues of B must all be positive for it to be 
guaranteed to be a completely positive map. If B has 
a negative eigenvalue but still transforms any positive 
p(0) into a positive p(t), then B is a positive but not 
necessarily completely positive map. 

Non-completely positive maps have been measured us- 
ing quantum process tomography (QPT) 105,106 which 
has caused the specifics of QPT to be questioned |107j . 
But the possibility that a map which is not a completely 
positive map can transform a valid quantum state into 
another valid state has brought a great deal of interest in 
studying the conditions for complete positivity. This is in 
addition to the interest in it due to the partial transpose 
as an indicator of entanglement pL 051 1109] , 

In 1994, Pechukas showed that complete positivity con- 
strains a system to product states of the form Pse = 
PS ® Pe, where pe is & fixed state of the bath [52 1 110] 
which excludes correlations and does not represent many 
physical situations. Alicki in Ref. [Ill] argued that there 
is no general definition for the reduced quantum dynam- 
ics beyond the weak coupling regime, therefore, when 
the system is in an initially correlated state with the en- 
vironment, linear assignment maps have no unique def- 
inition |107j . and linearity would only be preserved for 
states that are invariant under the transformation |lllj . 
Pechukas replied in Ref. |110j . and agreed that open 
system reduced dynamics can be non-linear. However, 
Rodriguez-Rosario et al. examine the assignment maps 
in |107j and argue against giving up linearity by noting 
that the assignment maps can be linear if the conditions 
of consistency or positivity are relaxed, and favor relaxing 
the positivity condition. A quantum system that inter- 
acts with the environment before our prescribed t = 
can be described by completely positive dynamics if the 
environment does not re-act on the system |104j . i.e. the 
coupling is weak and/or the initial state is in a particular 
form [9"I] , 



the Kraus representation or operator-sum decomposition 
|112j . although it was originally discussed in this con- 
text by Sudarshan, Mathews, and Rau [95] . Jordan, et 
al. demonstrated that entanglement in the initial state of 
the system can lead to non-completely positive maps that 
still transform a positive p into another positive p |113) . 
Rodriguez-Rosario, et al. found that for purely classical 
correlations, the "quantum discord" (defined below) van- 
ishes, and this is a sufficient condition for completely pos- 
itive reduced dynamics [114] . Later, Shabani and Lidar 
demonstrated that the quantum discord was also a nec- 
essary condition for complete positivity |115j . Quantum 
discord was introduced by Ollivier and Zurek in 2001, it 
is defined as a 'measure of the quantumness of the corre- 
lations' (116] . and is calculated as follows: 



S(S : E) =-Tr (p E log(p E )) + Tr (p SE log(p SE )) 



where H{x) — H(p x ) = — Tr (p x log(p^)) is the Von Neu- 
mann entropy and — J^j Tr(Hf pse) "/(n^pse) * s ^ e con ~ 
ditional entropy, defined as the entropy of the system 
with respect to a set of projective measurements per- 
formed on the environment. Quantum discord provides 
a measure of the nature of correlations, it vanishes for 
classical correlations and is maximum when there is en- 
tanglement. 

II. BACKGROUND 

As mentioned before, the extent to which the noise 
from the environment can be included in a quantum sim- 
ulation is dependent on both the simulating and simu- 
lated systems. Of course it would useful to have some pre- 
vious knowledge of the system-bath interactions. How- 
ever, this is often not the case. Here we study effects 
of unwanted noise in a quantum simulation using an 
algorithm [6] that simulates the one dimensional Fano- 
Anderson model. In this case we have a realistic model 
of the interaction and use the dynamical maps of the 
system to describe the noisy evolution. Starting with 
different initial states of the system and bath, we reduce 
the dynamics to a two-particle model system. The al- 
gorithm requires the two particles to be initialized in a 
particular state. Due to interactions with external qbits 
in the simulating device, these initial conditions may be 
imperfect. In addition, if the particles are allowed to 
interact for some small time before the begining of the 
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actual algorithm, the particles could begin in a correlated 
or entalged state. We consider the possibility of errors in 
the preparation of one of the particles in the system as 
well as the possibility of correlations between particles. 
We added a visualization of the evolution of the Bloch 
vector in order to provide an intuitive picture of the dif- 
ferences in the initial states and how they evolve. It is 
useful to note that, regardless of the non-complete posi- 
tivity of some of the maps obtained, the final state is a 
physical state and the system is a realistic physical model 
with realistic couplings. The significance of these results 
will be discussed in the conclusions. We now describe our 
methods and results. 



A. Quantum Algorithm 

Ortiz, et al. proposed an algorithm for the quantum 
simulation of the one-dimensional Fano- Anderson model 
[T5] , This model consists of an impurity described by 
an energy e surrounded by a ring of n spinless fermions 
having energies e ki . The fermions interact with the impu- 
rity, which is also a spinless fermion, through a hopping 
potential V [TS] . The diagonalized wave-number rep- 
resentation of the Fano-Anderson Hamiltonian is given 

by [61 ng 



to Eq. (TT8J) [B] 



xe-'i^e-^^e'^'^e-'i^e'^v. (20) 

The goal is to see if the initial state of the impurity has 
changed over time and, if so, how much. For this purpose, 
we use the time correlation function C(t) = b(t)b(0y , 
which in spin operator representation becomes C{t) — 



e iHt a l_ e -iHt a l + 



where a + = a x 



ioy and cr_ 



io y . The time correlation function provides information 
about the overlap of the initial and final states of the 
impurity. 



We use the same form of the Hamiltonian in Eq. ( 17 ) to 



perform the unitary evolution on different initial states of 
the system and perform the same operation regardless of 
prior interactions. We then obtain the reduced dynamics 
of the state of the impurity site (qbit 1) and then obtain 
the dynamical map that describes the evolution. We also 
calculate the time correlation function for the purpose of 
comparing the results of the different situations to those 
of an ideal scenario. In this way we observe the effects 
of the noise and possible errors in the outcome of the 
simulation. 



B. Simulation with Noise 



n-l 



H = Ye k J k c ki + etfb + VY(c{.b + tfc ki )6 ki0 . (16) 



i=0 



i=0 



The system is mapped via Jordan- Wigner transformation 
to the spin system to obtain [3] 



H 



£ k 2 , V 



. . 2 l "^ + a y°y)- ( 17 ) 



Ortiz, et al. consider an NMR device for their simula- 
tion as do we, but the model is certainly not limited to 
this type of device. The simulator has an NMR drift 
Hamiltonian of the form [3] 



(18) 

The control Hamiltonian for spins in the system is 

= X!^ * + "ft^l a i,J a l a i' ( 19 ) 




To include other qbits in the environment surrounding 
the system of interest we modified the control Hamilto- 
nian in two different ways: 

1. First, we added two spins and had them interacting 
via zz coupling with the particle that represents the 
state of the fermion site (qbit 2): 



NMR 



1 | (e + £fco) 



'e-£feo) 



V 2 



1 / (e + £fc ) , ( 0- £ &o) 



J zz 2 3 i "zz 2 4 i 3 4 

—£ a z a z + ~^- (T z (T z + -^^z^z- 



(21) 



2. Next, we added an extra particle, which interacts in 
the same fashion (zz coupling) with both particles 
that represent the system of interest: the resonant 
impurity and the fermion site: 



H 



NMR 



1 / ( £ + £feo) \( £ -£fc ) 



where the a are controllable. The last term is considered 
controllable because it can be turned on/off with the x 
and y rotations. 

To obtain the representation of the Hamiltonian in 



Eq. (17), the following control sequence can be applied 




6 



6, 



(a) t=0 



(b) t=0.3 





(c) t=0.6 



(d) t=0.9 



FIG. 1: Evolution of the Bloch Vector of the reduced 
dynamics of qbit 1 in the initial state p\ = |0) (0| as a 
function of time. 



where J zz represents the zz coupling constant. Wc 
used the same control sequence from Eq. (p0| to obtain 



Eq. (17), to represent a situation in which the extra qbits 



are environmental and thus are taken to be unknown. 



A. 2 Entangled states 

|V(0)> =a |01} + ai|10), (24) 

where aig + a\ — 1, and the density operator is 
given by p(0) = |V(0)> <V(0)|- 

A. 3 Correlated states 



P (o) = (i - p) (pi® ®pI 



(25) 



where p[ and p\ are the density operators cor- 
responding to some initial state of the impu- 
rity ( "spin-down" /occupied) and fermion ("spin- 
up" /unoccupied), respectively, and p] 1 and p 1 ^ 
correspond to the other initial state of the impu- 
rity ( "spin-up" /unoccupied) and fermion ( "spin- 
up" / unoccupied) . 

We represented the initial state of the impurity in terms 
of its x, y and z projections of the Bloch vector. The 
magnitude of each component of the projections, a%, can 
be obtained by performing the partial trace over every- 
thing else except qbit 1, as dj = Tr[o-j(p(0))]. 
First consider an initial density operator 



Ps(0)= *(1 



o5) 



In this case, case A.l, 



III. RESULTS 

In this section we describe the results of the simulations 
for the two different modifications to the Hamiltonian as 
well as different initial states. 



A. States with Bloch vector in the z direction 

We first consider states with only a z component to 
their Bloch vectors. These form a special class of states 
due to the commutativity of the zz Hamiltonian with 
these initial states. 



1. Noiseless Quantum Simulation 

Here we consider the cases where no bath is present, 
but different initial states are considered. Three cases 
are considered corresponding to three types of different 
initial states used in the simulation: 

A.l Pure states 

M0)> = |00},|01),|10), |11). (23) 
Density operator calculated as p(0) = \tp(0)) (V'(O)I- 



Ps(0) = -(l + a 3 a z ), 

where 03 represents a real constant that is equal to, or 
less than, the radius of the Bloch sphere (i.e. < 03 < 1). 
It represents the projection along the z axis. The final 
state was obtained through the reduced dynamics of ps 
after the evolution: 

ps(t)=Ti[p s (0) (Up(fl)U*)]. 

When the initial states ps(0) only had a z component, 
the final states ps(t) only had a z component as well 



Ps(t) = -(l + b 3 a z ), 

where 63 is another real constant that is subject to 
< 63 < 1. The value of 63 depends on a 3 and on 
the parameters e, e^, V and t. When states with only a 
z component are input, the final states also only have a 
z component. This is consistent with the hopping model 
where the "spin-down" corresponds to the state being 
occupied. The evolution is described by the dynamical 
map 



B 



( 


1+63 

2 








\ 







l+b 3 
2 
















1-63 

2 





V 











1-63 J 

2 / 



(26) 
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A. 6 Correlated states 




FIG. 2: Eigenvalues of the dynamical map of the 
reduced dynamics of qbit 1 in the initial state = 1 01) 
in the closed system. The parameters are e = 8, e = — 2, 
V = 4, for the time interval At € [0.1, 0.9]. 



The the eigenvalues of the map are plotted as functions 
of time in Figure ([2|. 

The case of a maximally entangled state case A. 2 also 
has a similar form with only a z component. Therefore 
the map is also given by a form similar to Eq. ( 26 1 



As 

can be seen in Fig. ([2]), the eigenvalues of the map cor- 
respond to a completely positive evolution, even though 
the quantum discord of the initial state was a maximum. 
Similarly, for case A. 3. We therefore note, for later ref- 
erence, that in these cases all states have only a z com- 
ponent in the initial and final states of the system. Thus 
there is only this standard interpretation of the hopping 
model Hamiltonian when there is no external noise. 



2. Simulation with noise from spin bath 

In this section we present the results for systems gov- 
erned by the Hamiltonians in Eqs. (21) and (22 1. The 



goal is to simulate a two body problem, so we used the 
same control sequence in Eq. (20 1. However, the initial 
state of a "bath" of two particles was included in the total 
system Hamiltonian. As in the simulation that had no 
external noise, we chose different initial configurations, 
Explicitly, including the bath qbits these are: 

A. 4 Pure states 

\t/>(0)) = |0011>, |0111), |1011), |1111), (27) 

and density operator p(0) = \tp(0)) (^(0)|. 

A. 5 Entangled states 

|-0(O)> =a |0111) +ai|1011) (28) 



Where q.q 



1, and the density operator is 



p(0) = ((1 - p)(p[ ® p{) + p{pY ® ^))®(|1) <1|)®(|1) (i|). 

(29) 

The fact that the states only had a component in the z 
direction and only interact with the bath via zz couplings 
resulted in results that were very similar to the ones in the 
previous section. The initial state of qbit 1 (the impurity) 
can again be written in Pauli notation as: 



ps(0) = Tr B p(0) = ^(l + a 3 a z ). 



(30) 



The final state is obtained by tracing over the bath de- 
grees of freedom 



Pi(t)=Tr E (Up(O)U^) 



(31) 



63 is another real constant that can be positive or neg- 
ative, depending on the direction of the Bloch vector of 
the final state along the z axis. 

The most general dynamical map has the same form 
as the map in Eq. ( 26 ) , 



B = 



1 


l+b 3 
2 
















1 + &3 

2 
















1-63 

2 





V 











l-b 3 
2 



\ 



(32) 



? iven by p(0) = |V>(0)> <V(0)| 



We observed that the coupling J zz has an effect in the 
rate of change of the state of qbit 1, which is presented 
in the results for the calculation of the time correlation 
function. For purposes of comparison, the parameters 
of the system were the same as the results above. In 
Figs. ([3]) and Q, the eigenvalues of B are plotted with 
the couplings to the particles of the spin bath being J zz — 
8 and J zz = jg respectively. 

Figs. ([2]), |3|) and Q show the evolution of the same 
initial state but each has a different environment. Being 
states initially in the z direction, the dynamics are com- 
pletely positive since the interaction with the bath is a zz 
coupling. However, it does change the hopping rate. In 
Fig. (|3| this is particularly noticeable due to the choice of 
the coupling. The state of the impurity does not transfer 
as easily due to the strong correlations generated by the 
interaction with the spin bath. In Fig. Q the situation 
is different. In this case the eigenvalues remained the 
same regardless of the strength of the coupling with the 
environment. Having a single extra particle interacting 
with both qbits with the same strength would mean they 
both have the same interaction with the bath. 



B. Arbitrary initial direction of the Bloch vector 

Noise in the initialization of the state could result in 
a direction for the Bloch vector which is not in the z 
direction. States that have an x or a y component to their 




FIG. 3: Eigenvalues of the dynamical map of the 
reduced dynamics of qbit 1 . Two additional qbits are 
interacting via zz coupling with qbit 2, J zz — 8. The 
initial state of the total system and bath is 

\ip) = |01 1 1) . The system parameters are e = —8, 
e = — 2 and V = 4, in the time interval t € [0.1,0.9]. 



Larmor precession will be present. This is often observed 
in a NMR device under general circumstances and leads 
to noise in the system. Here we consider an initial state 
with a component of the Bloch vector in the x direction. 
Clearly a y component is not necessary, and only specifies 
a different initial condition for the angle since the system 
will precess. 



1. Noiseless Quantum Simulation 

The initial states were chosen to have a component in 
the x direction; the components in x and z were selected 
such that the magnitude of the Bloch vector is close to 1 
emulating a small error in the initialization. Explicitly, 
the different initial configurations were: 

1 



Pl(0) = o(l + ai<Jx + a 3&z), 



P2(0) = -(1 + a\o x - a 3 a z ), 




FIG. 4: Eigenvalues of the dynamical map of the 
reduced dynamics of qbit 1. An additional qbit is 
interacting via zz coupling with qbits 1 and 2, 
J zz = 1/10. The initial state is |V>) = |011). The system 
parameters are e = —8, e = —2 and V — 4, for times 

t e [0.1,0.91. 



polarization vector, or Bloch vector, exhibit precession 
and approximate more accurately what happens in a real 
experimental situation. In this case, if some initial error 
on the system has a component in the x or y direction, 



where pi is the state of the impurity, p2 is the sta te of 
the fermion and the aj are subject to < \/a\+ aj < 1. 
Therefore, the total initial state is 

/9(0) = pi(0)®p 2 (0). 

The final state of the impurity was, once again, ob- 
tained by doing a partial trace over the degrees of free- 
dom of the fermion 

p{t) = Tr B (Up(0)tf) - ^(l + hax + hay + ha,). (33) 



The map B is given by 



B 







/ 

Q 6l_ 1-^3 

Oi 2 



-ib 2 \ 

o 





2 



(34) 



/ 



The eigenvalues of B are given by 



Ai = 

A 3 = 
where 



^Abj+albj 



2a x 



2oi 



, A4 — 



oi + ^Ab\ + a\bl 

20! 



a 1 



2ai 



(35) 
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b t = { cos ( -t(e + e ko )j cos (±ty/4V 2 + (e-e ko )A - sin (~t(e + e ko ) 



b 2 ={- sin ( h(e + e ko )\ cos fyy/4V 2 + (e - e ko ) 2 J - cos ( ~t(e + e feo ) 



(e - e) sin (§ V^ 2 + (e - e feo ) 2 ) 
( 6 _ e ) S i n (i t ^ + ( € - £fco )2) 



ai, 



and 



6a = 



J(-l + a 3 )^ 2 + a 3 (e - e? + (1 + a 3 )^ 2 cos (\t^W 2 + (e - e) 2 ) 

4F 2 + (e - e) 2 



Oi 



(36) 



Note that if ai >-> 0, then &i and 62 are 0. The factor 
a\ in the denominator of the eigenvalues is eliminated 
using l'Hospital's rule, and that yields 

I-63 . I + 63 



Ai = 
A3 = 



2 

I-63 



A 2 = 
A 4 = 



1 

I + 63 



(37) 



which are the same as the eigenvalues of the map in 
Eq. ( 26 ) . The eigenvalues of B when ai > are shown 
in Figure In Fig. ([6]), the dynamics of the system are 
positive but not completely positive. This system is not 
in contact with a bath or reservoir, but it consists of two 
particles. This is a case of initial correlations between 
particles in the system, which are errors for this model 
since correlations should not be present in initial state 
preparation. The general observation was that when the 
initial state has a component of the Bloch vector in x 
or y as well as one in z, the result is a non completely 
positive map. 



2. Simulation with noise from spin bath 

The results in this subsection are generated from 
adding the qbits in the spin bath, and using the following 
initial states 



p(0)=/»i(0)®P2(0)®(|l)(l|)®(|l)<l|), 



where 



and 



Pl(°) = o(l + + «30- z ), 



Pz(0) = ^(1 + ai<r x - a 3 a z ). 



(38) 



(39) 



(40) 



The reduced dynamics of S are given by 



1 



with a B map of the same for as that in Eq. (34), 

\ 



B = 











— ib 2 
ai 





l+b 3 


61 








2 

Zll 


2 





\ ai 


01 






1-63 

2 



(42) 



p(t) = Tie (U P {0)rf) = -(l + b 1 a x + b 2 a y + b 3 a z ), (41) 



Once again, the noise, which has the form of purely zz 
couplings, caused variations in the parameters, mostly in 
the rate of change of the state of qbit 1. The eigenvalues 
for a system with two spins interacting with the fermion 
only and for one spin interacting with both particles in 
the system are presented in Figs. ^ and 

In Figs. ([7]) and rtSJ) the reduced dynamics are not com- 
pletely positive. This is due to the initial state of the 
impurity site (qbit 1) having a component of its Bloch 
vector in the x direction. The algorithm was designed 
to have an initial state where one of the two state sys- 
tems is in the up state and the rest are in the down 
state. Dynamical maps obtained through quantum pro- 
cess tomography can present discrepancies if the initial 
states are prepared through different experimental meth- 
ods [90j . Thus the x component represents a preparation 
error which gives rise to a non-completely positive map 
like in the previous case. 



C. Time correlation function 

Ortiz et al. calculated the time correlation function 
C(t) = 6(t)6(Q)t, and plotted the result as |G| 2 = 
Tr (p(t)p(O)) as a function of time. Since we want to cal- 
culate the effects of noise and different initial state, we 
followed the same procedure for the different situations. 
The results are summarized in graphs, Figs. Q, (10) 
and ( |TT| ). In Fig. there is a slight difference between 
the results of the original system compared to those un- 
der which errors could arise due to noise and unknown 
initial states. The coupling to the environment affects 
how fast or slow qbit 1 evolves. However, if the coupling 
to the bath is weak, these errors are not as prominent. 
There was one case in which there was no effect on the 
speed of change by spins in the bath. 

When the initial state had a component in x the result- 
ing correlation functions were very close to the original 
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(a) t=0 




(c) t=0.2 



6, 



(b) t=0.1 




(d) t=0.3 




FIG. 6: Eigenvalues for the dynamical map of the 
reduced dynamics of qbit 1. The initial state of the 
system is p x = \ (1 + 0.2a x + 0.97a z ), p 2 = \ (1 - a z ). 
The parameters of the system are e = —8, e = —2, 
V = 4. Evaluated in the time interval t € [0.1,0.9]. 



(e) t=0.4 



v / 



(f) t=0.5 




(g) t=0.6 



(h) t=0.7 



FIG. 5: Animation of the evolution of the Bloch Vector 
of the reduced dynamics of qbit 1 in the initial state 
pi = 1(1 + 0.2^ + 0.97^) 



problem. This is important because a small error like 
this one may not be easily identified in the time corre- 
lation function. In Fig. ( 10 ) we show how the coupling 



to a spin bath can affect the rate of change of the evolu- 
tion. As mentioned before, these results only include zz 
couplings. The strength of the couplings were adjusted 
in order to see the effects. 

Because quantum simulations are performed on quan- 
tum systems, where access to complete information about 
the state at all times is not available, correlations with 
the bath can be by detected by differences in the rate 
of change of the evolution. However, it is useful to also 
study the case in which the coupling is not only in z. In 
Fig. (11), we increased a x , the component of the Bloch 
vector in x, to see how it affects the final result. When 



FIG. 7: Eigenvalues of the dynamical map for the 
reduced dynamics of qbit 1. The initial state of the 
system is p\ = | (1 - 

P3 = |(l-^) 

parameters are e = —8, e = —2, V = 4, J zz — ^ 
evaluated in the time interval t £ [0.1,0.9]. 



0.2^ + 0.970, P2 = |(l-o-,), 
Pi = | (1 — <J Z ). The system 



the x component of the Bloch vector is increased, we can 
see shifts in the time correlation function. The greater 
a% is, the larger the observed shift. This could be useful 
for detecting possible errors in state preparation. 



IV. CONCLUSIONS 

Interactions of quantum systems with a surrounding 
environment are undesirable for reliable quantum simu- 
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FIG. 8: Eigenvalues of the dynamical map for the 
reduced dynamics of qbit 1 . The initial state of the 
system is pi = \ (1 + 0.2a x + 0.97cr z ), p 2 = \ (1 - o z ) , 

p3 = o (1 — Pz)- The system parameters are e = — 8, 
e = —2, V = 4, J zz = evaluated in the time interval 
t £ [0.1,0.9]. 




0.2 0.4 0.6 0.S I 

Time (a.u.) 



FIG. 10: Time correlation function of the reduced 
dynamics of qbit 1. The system parameters are e = —8, 

e = -2, V = 4 in the time interval t € [0.1, 0.9]. The 
results correspond to the closed system and the system 
that interacts with two additional qbits, coupled only to 
qbit 2. The initial state ofqbit 1 is p = |0) (0| for one set 
of results, and p = | (1 + 0.2a x + 0.97a z ) for the other. 




FIG. 9: Time correlation function of the reduced 
dynamics of qbit 1. The system parameters are e = —8, 
e = -2, V = 4 and time interval t G [0.1,0.9]. The 
results represent the closed system, the system where 
qbit 2 interacts with two additional qbits, the system in 
which an additional qbit that interacts with qbits 1 and 
2. This was done when qbit 1 was in the initial states 
p = |0) (0| and p = |(1 + 0.2a x + 0.97cr z ), as indicated 
above. 



FIG. 11: Time correlation function of the reduced 
dynamics of qbit 1. The system parameters are e = —8, 

e = —2, V = 4 evaluated in the time interval 
t £ [0.1,0.9]. The result represents the time correlation 

function of the closed system compared to the 
correlation function of the reduced dynamics of qbit 1 
in the initial state p = ^ (1 + a\o x + a^a z ) for different 
values of a\ and a^. 



lations and for quantum information processing in gen- 
eral. Understanding and controlling or suppressing the 
noise from the environment is one of the most important 
objectives of studying open system quantum dynamics. 
Lloyd's suggestion to use the noise to simulate the in- 
teraction of the system with the environment is clearly 



useful only in special cases. For some analog simulators, 
isolation has nearly been achieved [IS], but that will not 
be the case for many devices. For the cases where noise 
suppression is required, understanding the noise will be 
necessary for implementing the appropriate noise control 
method. 
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It is known that interactions with the environment can 
lead to correlations that can result in non completely pos- 
itive maps. We found that such non completely positive 
maps are not rare in our study of a very simple model of 
a quantum system of fermions which can readily be sim- 
ulated on a quantum computing device, or a dedicated 
quantum simulator. This Fano- Anderson model exhibits 
maps which are not completely positive for a variety of 
initial states, some of which were entangled and some 
with other non-trivial quantum correlations in the sense 
of non-zero quantum discord. They were shown to arise 
for even a fairly small transversal component to an ini- 
tial density matrix which is supposed to have its Bloch 
vector aligned along the z axis. Thus fairly small experi- 
mental errors can lead to maps which are not completely 
positive in a rather simple experiment. These noises also 
cause relatively small errors in the final outcome of the 
measurement. 

Initially correlated states, if they are not so identified, 



but are instead identified improperly as arising from com- 
pletely positive maps, may encourage an experimenter to 
try to employ dynamical decoupling controls to eliminate 
errors. These controls will be ineffective in these cases. 

We have used a very specific and simple model to illus- 
trate the effects of noise on the system including the pres- 
ences of maps which are not completely positive. How- 
ever, it is important to emphasize that these effects are 
quite general and will be present in some form in many 
other quantum systems including a wide class of quantum 
simulations. 
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